Power-Law Behavior of Bond Energy Correlators in a Kitaev-type Model with a Stable Parton 

Fermi Surface 



O 

(N 

Oh: 

< 



I 

o 
o 



> 
in 



Hsin-Hua Lai and Olexei I. Motrunich 

Department of Physics, California Institute of Technology, Pasadena, California 91125, USA 

(Dated: April 14, 2011) 

We study bond energy correlation functions in an exactly solvable quantum spin model of Kitaev type on the 
kagome lattice with stable Fermi surface of partons proposed recently by Chua et al., Ref. | arXiv: 1010. 1035 1. 
Even though any spin correlations are ultra-short ranged, we find that the bond energy correlations have power 
law behavior with a 1/1^^ envelope and oscillations at incommensurate wavevectors. We determine the corre- 
sponding singular surfaces in momentum space, which provide a gauge-invariant characterization of this gapless 
spin liquid. 



I. INTRODUCTION 

In the last two decades, there has been dramatic theoret- 
ical progress in our understanding of RVB ideas' and spin 
liquids. We now know that there are many different kinds 
of spin liquids. Gapped topological spin liquids^— are best 
understood and have been shown to exist in model systems. 
Gapless spin liquids are also possibl o^-'"''^ and recently real- 
ized in experiments^^— , but are understood to a lesser degree, 
particularly when both the emergent parton and gauge field 
excitations are gaplessi^i^l^ 

From the early days, slave particle approaches^ have played 
an important role in studying such phases. The discovery by 
Kitaev '^^ of an exactly solvable, interacting two-dimensional 
spin- 1/2 model on the honeycomb lattice with a spin liq- 
uid phase paved an exciting road for the study of spin liq- 
uids. Since then, there have been many studies of Kitaev-type 
modelsi^— . In the last year, several of them realized gapless 
spin liquid phases with parton Fermi surface^/^^^^^ (and 
gapped Z2 gauge fields). In this work, we want to directly de- 
tect the presence of a surface of low-energy excitations. Note 
that the parton Fermi surface itself is gauge-dependent and is 
not accessible via local observables. However, there is a geo- 
metric surface information that is physical and can be detected 
using gauge-invariant local energy observables. We take up a 
very recent model by Chua et ah^^^ to illustrate this point. 

Chua et al^ proposed an exactly solvable spin-3/2 model 
on the kagome lattice and found a regime with a gapless spin 
liquid with a stable Fermi surface. Motivated by such a spin 
liquid phase and known techniques to characterize situations 



with gapless partons 
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we propose to study gauge- 



invariant operators such as bond energy operators. Our main 
results show that, unlike any spin correlations which are ul- 
tra short-ranged, the bond energy correlations have power-law 
behaviors with 1/|t"|'^ envelope in real space and also oscil- 
lations at incommensurate wavevectors which form what we 
call singular surfaces^^^^^ in the momentum space. An in- 
teresting aspect of the Chua et al^ model is that there is no 
"nesting" in the Cooper channel for the low-energy fermions 
because of the absence of inversion symmetry and the broken 
time -reversal. This gives a non-trivial fe^fi + ^fl critical sur- 
face in addition to more familiar kpR — kpL (a.k.a. "2kp") 
surface in the local energy correlations. 

In connection to experiments, the physics discussed here 



can be conceptually related to the recent gapless spin 
liquids in the organic compounds K-(ET)2Cu2(CN)3and 
EtMe3Sb[Pd(dmit)2]2r^"— where much thinking focused 
around the possibility of gapless Fermi surface of spinons. Al- 
though the Kitaev-type theoretical models are not directly ap- 
propriate for these materials, some of the qualitative physics 
discussed here applies more generally to gapless spin liq- 
uids and has practical implications. Thus, such Ikp physics 
information can also be revealed by measuring Ruderman- 
Kittel-Kasuya-Yosida (RKKY) interaction between magnetic 
impurities^2iS or by measuring textures in local spin suscepti- 
bility (Knight shift experiments) or other local properties near 
a non-magnetic impurity.i^^i^ We will discuss further possible 
connections in the conclusion section. 

We also mention that entanglement properties of a ground 
state wavefunction can be used for characterizing a phase 
of matter, especially for gapless spin liquids, in addition to 
gauge-invariant observables with power law correlations. For 
instance, a recent paper^- measured entanglement entropy in 
the Gutzwiller-projected Fermi sea wavefunction on the trian- 
gular lattice and found logarithmic violation of the area law, 
which strongly suggests the existence of gapless Fermi surface 
in the resulting spin liquid state. Other recent works^ used 
the entanglement entropy to estimate the central charge in 
DMRG studies of spin- 1/2 Hamiltonians with ring exchanges 
on multi-leg ladders, and found that the central charge in- 
creases with the number of legs as expected in such gapless 
spin liquids. 

The paper is organized as follows. In Sec. we start 
from the Chua et al. HamiltonianSS, on the kagome lattice. 
In Sec. Hn] we define bond energy correlation function. In 
Sec, nil Al we provide a theoretical approach to describe the 
long-distance behavior of the coiTelations. In Sec. HUB] we 
present exact numerical calculations of the bond energy corre- 
lations. We conclude with some speculations about similarity 
with recent experiments in EtMe3Sb[Pd(dmit)2]2 in which a 
gapless spin liquid has been realized J^— 



II. CHUA-YAO-FIETE KITAEV-TYPE HAMILTONIAN 

We begin by formulating the Hamiltonian in the parton 
language. The model is defined on the kagome lattice, see 
Fig. [T] On each site i of the kagome lattice, there is a phys- 
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ical four-dimensional Hilbert space realized using six Majo- 
rana fermions (^f, (^f, Ci, and di, with the constraint 
Di = —i^l^fQ£,fcidi = 1 (namely, for any physical state 
|$)phys, we require A|'J')phys = |^')phys)- The Chua-Yao- 
Fiete Kitaev-type Hamiltonian is 

■H = i ^ Uij [JijCiCj + J[jdidj\ + iJa ^ Cid^ (1) 

(y) » 

-a^ VKo - /? J2 W^aW^v , (2) 

O (A,V> 

where (ij) represents nearest neighbor links and Uij = —Uji, 



with Uij = if {ij) G A and ; 



for bond directions chosen to go counter-clockwise around the 
triangles. Placket operators Wp = Y[{ij)ep'^'^j' with p = 
A, V, O, are gauge-invariant (i.e., act in the physical Hilbert 
space) and are conserved by the Hamiltonian. The terms in 
Eq. (|2]i with a > and /3 > are added to stabilize particular 
ground states with Wo = 1, Wa = = ±1- Since in the 
Kitaev-type model, [uij ,7{\ = [ ] = 0, we can treat 

the Z2 gauge fields as static background and replace by 
their eigenvalues ±1. We then have free Majorana fermions 
c and d hopping on the lattice in the presence of "fluxes" (j)p 



defined via e 



n 



Throughout, we work in the ground state with Wq = 1, 
Wa = W\7 = 1 which breaks time reversal symmetry; this 
translates to fluxes {0o, '/'A, '/'v} = {^r, 7r/2, n/2} as shown 
in Fig. [U We fix the gauge by taking Uij ~ 1 with bonds 
i ^ j directed counter-clockwise around the triangles. There 
are three physical sites per unit cell and six remaining Ma- 
joranas per unit cell. We replace the labeling {ci,di} with 
^' r 1 , where r runs over the Bravais lattice of unit cells 
of the kagome network and a runs over the six Majoranas in 
each unit cell (three c Majoranas and three d Majoranas). The 
Hamiltonian can be written in a concise form. 



{{r,a),{r',a')) 



^r.a^r,a; r' ,a' ^ r'a' 



(7,7) 



(3) 



(4) 



The Majorana field satisfies the usual anticommutation rela- 
tion, = 25rr'Saa' = 2J/7- Inthcchoscn 
gauge, there is translational symmetry between different unit 
cells; hence, Ar,a- r'y = Aaa'{r - r'). 

In order to give a concise long wavelength description, it 
will be convenient to use familiar complex fermion fields. To 
this end, we can proceed as follows. For a general Majo- 
rana problem specified by matrix Ajj, we diagonalize Ajj 
for spectra, but only half of the bands are needed while the 
rest of the bands can be obtained by a specific relation and are 
redundant. Explicitly, for a system with 2m bands, we can 
divide them into two groups. The first group contains bands 
from 1 to in with eigenvector-eigenenergy pairs {vt.k, ^b.k}, 
where = 1,2, ... ,m are band indices, and the second group 
contains bands from m + 1 to 2m related to the first group, 

{Vb'=m+b.k, ^b' = m+b,k} = ~£6,-fc}- 

In the present case, 2m ^ 6 and therefore three bands are 
sufficient to give us a full solution of the Majorana problem. 



2^ /I 2^ TT /V 2^ '1 2^ /1 



FIG. 1 . Kagome lattice with three sites (labelled as 1,2, and 3) per 
unit cell. We consider Chua et al^ model with the ground state flux 
configuration as shown, {0o,<^a,<^v} = {^i", 'r/2, 7r/2}. We fix 
the gauge by taking iuij — i with bonds directed counter-clockwise 
around the triangles. We also show the bond energy operator Si2c 
whose correlations are presented in this paper, while other bond en- 
ergy operators have qualitatively similar correlations. 



For an illustration of how all 2m = 6 bands vary with mo- 
mentum k, we show the six bands in Fig. [2] along a cut with 
ky = 0. We label the bands from top to bottom as 1 to 6, 
and only bands 1 to 3 are used for the solution of the Majo- 
rana problem. Specifically, we write the original Majoranas in 
terms of usual complex fermions as 



*/' = E E V2[vbMi)Mk)+v,:^ii)fiik) 

b=l keB.Z. 
^ 3 



= \/]^E E [e*Xfc(«)/fc(fc)+H.c.] , (5) 



b=l keB.Z. 



where we used Vt^fc(/ = {r,a}) 



the number of unit cells), and the complex fermion field / sat- 
isfies the usual anti-commutaion relation, {/{[(fe), /fc(fe')} = 
Sbb' Skk' ■ In terms of the complex fermion fields, the Hamilto- 
nian becomes 



3 

^ = E E et,k\2fi{k)Mk)-i 

b=l keB.Z 



(6) 



Considering the model parameters from Chua et al^ real- 
izing the spin liquid with stable Fermi sea, we find that among 
these three bands, only ^ crosses the zero energy. Hence, 
as far as the long-distance properties are concerned, we can 
retain only band 3 and its Fermi surface is shown in Fig.|3] 



III. BOND ENERGY CORRELATORS 

In this paper, we focus on the bond energy correlation func- 
tions. There are several distinct bond energy operators one 
can consider However, all of them have similar long-distance 
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FIG. 2. Illustration of the energy spectra of the 2m — 6 bands 
along a cut with ky — 0. Here we take the same parameters as 
in Chuaefa/.j££ {Ja, Jv, Jk, J'v. J's} = {1.0,0.3,0.8,0.5,1.4}. 
Solving the Aij matrix in Eq. Q, there are 2m = 6 bands which we 
label from top to bottom as I to 6. Only half of the bands - e.g., I, 2, 
and 3 - are used to solve the Majorana problem, as long as the others 
- 4, 5, and 6 - can be obtained via the relation between eigenvector- 
eigenenergy pairs such as {v4,k,'^4.k} = {«3,-fc> -es.-fe}- 




Geir) = (Si2c(0)Si2c(r)) - (Bi2c(0)>(Si2c(r)) . (8) 

Power-law correlations in real space correspond to singulari- 
ties in momentum space, which we can study by considering 
the structure factor 



SB{q)=Y.GB{r)e 



(9) 



We will present exact numerical calculation of the bond en- 
ergy correlations in Sec. lIII Bl using the definitions in Eqs. O- 
Before showing the numerical data, we present a long 
wavelength analysis of such correlations due to the gapless 
Fermi sea of partons. 



A. Long wavelength analysis 

Focusing on the long distance behavior and therefore re- 
taining only the contribution from band-3, the bond operator, 
Eq. (ITJ, can be written approximately as 

6i2c(r) ^ J2 { [^^fefe'/3(fc)/3(fc')e'^'=+'='^ '' + H.c 

fe,fe'6B.Z. 

+ [Nkk'fl{k)f3{k')e 



-i{h—h')-r 



H.c. 



,(10) 



where Mkk' = 2i JAU3.fc(l)w3,fc' (2)/iVu 
2*Awlfc(l)w3,fc'(2)/iVuc. 

In order to determine long-distance behavior at separation 
r, we focus on patches near the Fermi surface of band 3 where 
the group velocity is parallel or antiparallel to the observation 
direction n = r/lrj, because at long distance \r\ ^ fc^^, the 
main contributions to the bond energy correlations come pre- 
cisely from such patches. Specifically, we introduce Right(R) 
and Left(L) Fermi patch fields and the corresponding energies 



fcfc' 



/<r)(<5fc) = /3(fe<r; + <5fc). 



~'fp\ 



(fl) 



PSh 



(11) 
(12) 



FIG. 3. Illustration of the Fermi pocket for true complex fermion 
fs in Eqs. (|5),l|6). The parameters are the same as in Fig. |2] The 
hexagon is the Brillouin zone boundary. Note that there is only one 
Fermi pocket for this set of parameters and the three pockets shown 
are related by a reciprocal lattice vector and are thus equivalent. 



behavior, so we present correlations for the bond energy op- 
erator corresponding to Ja term in the Hamiltonian between 
sites 1 and 2 as indicated in Fig.[T|and defined as 

Bl2c{r) = iJAUijCtCj = iJACr.lCr,2 = iJA*ilfi*^f2,(7) 

where from the first to second equation we specified to the 
working gauge. We will study bond energy correlator defined 



where the superscript (n) refers to the observation direction 

and P = R/L = +/— ; v^^p is the corresponding group ve- 
locity (parallel to n for the Right patch and anti-parallel for 
the Left patch); Q!p=/?/L is the curvature of the Fermi surface 
at the Right/Left patch; Jfey and dk± are respectively compo- 
nents of 5k parallel and perpendicular to n. It is convenient 
to define slowly varying fields in real space 



E 



(13) 



(5fcG Fermi Patch 



which vary slowly on the scale of the lattice spacing [and from 
now on we will drop the superscript (n)]. Therefore, in this 
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long wavelength analysis, the relevant terms in the bond oper- 
ator are 



BMr) - [{Nrr + N*^^j,) /I(r)/fl(r) + {R ^ L) 
+ [{Nlr + N*j,l) .fl{r)fR{r)e^^^---^--^-^ + H.c. 
+ \{Mrl - Mlr) .fe(r-)/L(r)e*('=-"+'=--)"- + H.c. 



(14) 
(15) 
, (16) 



where we dropped terms such as fR{r)fR{r) due to Pauli 
exclusion principle. The above long wavelength expression 
for the bond energy operator implies that the corresponding 
correlation function defined in Eq. ([8]l contains contributions 
with g = 0, ±{kpR — kp^), and ±{kpR + kp^). 

More explicitly, for a patch specified by ep{Sk) above, 
Eqs. (fTni-(fT2]i. we can derive the Green's function for the con- 
tinuum complex fermion fields as 



{fl./LWfR/L{r)) 



exp[Ti^] 



23/2^3/2^V2^|^|3/2 



(17) 



Using this and the long-wavelength expression for bond en- 
ergy operators, we can obtain the bond energy correlation, 
Eq. (IDl, 



GBir) 



(Nrr + N*^^^)^ [NLL + Nl^f 



aR\r\ 



aL\r\ 



2\N, 



RL 



nIr^ 



iiin[{kpR - kpp) ■ r] 



1/2 1/2, |Q 



2\Mrl - Mifl|2 cos[{kFR + kpp) ■ r] 



1/2 1/2, 



(18) 
(19) 

(20) 



Therefore, the above low energy description can be used 
to analyze the numerical data we obtain by exact calculations. 
Here we also note that the model does not have inversion sym- 
metry (and the time reversal is broken in the ground state), 
so the location of the corresponding R-L patches which are 
parallel or antiparallel to the observation direction can not be 
determined easily and need to be found numerically. 



B. Exact numerical calculation 

We calculate the bond energy correlations, Eq. (|8]l, for any 
real-space separation r and confirm that they have power law 
envelope l/lrp. For an illustration, we show the bond en- 
ergy correlations for r along a specific direction, e.g. x-axis, 
calculated on a 300 x 300 lattice. In Fig. |4] the log-log plot 
of |GB(r)| along the x-axis clearly shows the l/|Tp enve- 
lope. In addition, the irregular behavior of the data is due to 
oscillating components. For certain directions, the oscillating 
parts are sufficiently strong that GB(r) also changes signs. 
The wavevectors of the real-space oscillations form some sin- 
gular surfaces in the momentum space, which we will analyze 
next. 

Shifting our focus on the structure factor Ssiq) defined in 
Eq. (|9|l, we calculate the bond energy correlation at each site 



l<fil2c(0)Bl2c(x)>| 



10- 



10- 



10- 



IQ- 




10 20 50 100 



FIG. 4. Illustration of power-law behavior of bond energy correla- 
tion function Eq. We calculate Gs(r) with r taken along the 
i-axis for a system containing 300 x 300 unit cells. The log-log plot 
clearly shows envelope (straight line in the figure). Here we 
show the absolute values |Gs(r)| and indicate the sign with open 
square boxes for negative correlations and filled circles for positive 
correlations. The irregular behavior is due to oscillating parts; these 
appear to be rather weak, but if we change the observation direction, 
the oscillating parts can be stronger. 



within a 100 x 100 lattice and numerically take Fourier trans- 
formation. Figure [5(a)] gives a three-dimensional (3D) view of 
the structure factor. We can clearly see cone-shaped singular- 
ity at <7 = 0, which is expected from Eq. ( fTSl l, 



SB{q-0)^\q\ 



(21) 



A closer look at the structure factor also reveals singular 
surfaces at kpp — kpp and kpp + kpp, as expected from 
Eqs. (fT9] l and ( |20] |. In order to see the location of the singu- 
lar surfaces more clearly and compare it with our long wave- 
length analysis, we show top view of Ssiq) in Fig. |5(b)| We 
also numerically calculate Q± = kpp zLkpp (by first finding 
corresponding Right and Left Fermi points with anti-parallel 
group velocities) and superpose these lines on the figure. We 
can see that the lines we get from the long wavelength anal- 
ysis match the singular features in the exact structure factor. 
Note that the singularities are expected to be one-sided, 



SB{Q-+Sq)^\Sq^^\'^^e{-6q^^), 
SBiQ+ + 6q) ^ \Sqn f^^O [-5qn sign(a^ 



(22) 
(23) 



The first line is singular from the inner side of the "ring" in 
Fig. |5(b)| and the second fine from the inner side of the "trian- 
gles". 



IV. CONCLUSION 

We studied bond energy correlation functions in the 
Chua et al^ Kitaev-type model with a parton Fermi surface. 
Unlike spin correlations, we found that the local energy cor- 
relations have power-law behavior in real space with an enve- 
lope of l/lrj'^ and oscillations at incommensurate wavevec- 
tors that form singular surfaces in momentum space. By 
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2 H 




4 n 
3 

(a)Thiee-dimensional view of the stmcture factor 
2 JT I 




-27T ±JL 2jt 

3 3 

(b)Top view of tlie structure factor 

FIG. 5. (a) 3D view of the structure factor of the bond energy corre- 
lation, 5s (<j), defined in Eq. ((9). We can clearly see the singularity 
Ss(q) ~ |q| at q = and we also see weak singular lines, one 
forming a closed ring, and additional lines near the corners of the 
Brollouin zone, (b) These singular lines are brought out more clearly 
when the structure factor is viewed from top. We also superposed 
the locations of the singularities calculated using the Fermi surface 
information: The inner red ring specifies the line at kpR — kpL and 
the outer blue triangles specify kpR + kpL lines. 



combining low-energy theoretical analysis and exact numer- 
ical calculations, we determined the locations of the singular 
surfaces. These bond energy correlations provide a gauge- 
invariant characterization of such gapless spin liquid. 

We conclude by speculating about some interesting simi- 
larity with recent experiments in EtMe3Sb[Pd(dmit)2]2— — 
While the thermal conductivity measurements^ are consis- 
tent with the presence of a Fermi surface of fermionic exci- 
tations down to the lowest temperatures, very recent NMR 
experiments^ show a drastic reduction in spin relaxation be- 
low temperature of the order 1 K, almost as if a spin gap is 
opened. This reminds of the present situation where the spin 
operators have short-range correlations, which occurs because 
some of the constituent partons have a gap (here are ultra- 
localized), while there remain partons that are metallic and 
give rise to metal-like thermodynamics and manifestly gapless 
properties such as the discussed local energy correlations. Of 
course, the present model is on a different lattice and is very 
differently motivated. However, in a recent paper^ - working in 
a setting closer to the EtMe3Sb[Pd(dmit)2]2 experiments, we 
discussed the following scenario in magnetic Zeeman field: 
Upon writing the spin operator as 5*+ = f^fi, we considered 
a state where one spinon species (say, /^) becomes gapped 
due to pairing, while the other species retains the Fermi sur- 
face. In this case, spin correlations are short-range while 
the thermodynamics is metal-like. Furthermore, just as in the 
present paper, there are other properties that are manifestly 
gapless, e.g., spin correlations and transverse spin-2 cor- 
relations. It would be interesting to explore such scenarios in 
more realistic settings further. 
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